Spin/Orbital Pattern-Dependent Polaron Absorption in Ndi x Sr x Mn03 
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We investigated optical properties of Ndi-^Sr^MnOg (x= 0.40, 0.50, 0.55, and 0.65) single crys- 
tals. In the spin/orbital disordered state, their conductivity spectra look quite similar, and the 
strength of the mid-infrared absorption peak is proportional to x(l—x) consistent with the polaron 
picture. As temperature lowers, the Ndi-^Sr^MnOa samples enter into various spin/orbital ordered 
states, whose optical responses are quite different. These optical responses can be explained by the 
spin/orbital ordering pattern-dependent polaron hopping. 

PACS number; 75.30.-m, 78.30.-j, 72.80.Ga, 71.38,+i 



Recently, lots of attention has been paid to the roles 
of orbital degree of freedom in the physical properties 
of manganites. Resonant x-ray scattering of a layered 
manganite provided evidences of orbital ordering [[l] . Po- 
larization microscope images demonstrated the existence 
of optical anisotropy , which might be due to the or- 
bital ordering. A new elementary excitation caused by 
symmetry breaking of orbital ordering, called 'orbiton', 
was proposed and claimed to be observed in LaMnC>3 by 
Raman spectroscopy measurements Since the cor- 
responding Raman peaks were observed in the energy 
region of a few tenth of eV, similar excitations could be 
observed by infrared spectroscopy. 

Indeed, an anomalous incoherent absorption peak has 
been observed in the mid-infrared (IR) region for nu- 
merous manganites (|-[|. However, the origin of this 
mid-IR absorption peak is not fully understood yet. In 
one scenario, it was suggested that the mid-IR absorption 
should be due to orbital excitations . Even though the 
observed anisotropy might be explained by this model, 
there have been few experimental reports to compare 
with the theoretical predictions even on a qualitative 
level ||. In the other scenario, it has been suggested 
that the mid-IR absorption should arise from hopping of 
the Jahn-Teller (JT) polaron [@-|Jll]Jl|l . The polaron 



binding energy obtained from the mid-IR peak position 
agrees with that from other measurements. However, this 
scenario has not been extended to explain the observed 
optical anisotropy. Till now, there is no consensus on 
these two scenarios. Therefore, it is important to un- 
derstand the role of the orbital degree of freedom on the 
mid-IR absorption peaks. 

Ndi_ K Sr x Mn0 3 (NSMO) shows various types of 
spin/orbital (S/O) ordering patterns at low temperature 
(T), depending on the value of the doping concentra- 
tion (x) |L3]]. Optical investigations on NSMO in various 
S/O ordered states can provide us further insights on the 



role of S/O for the mid-IR absorption in manganites. In 
spite of such importance, few investigations have been 
reported due to the difficulties in growing high quality 
single-domain crystals. 

In this Letter, we present optical responses of NSMO 
( 2=0. 40, 0.50, 0.55, and 0.65) multi-domain single crys- 
tals and a single-domain crystal with 2=0.55. At room 
T, optical conductivity spectra cr(o>) of all the samples 
appear quite similar and can be explained within the 
JT polaron picture. At low T, by taking account of 
the effects of the S/O ordering pattern on the polaron 
hopping motion, the anisotropy in cr(u>) for the single- 
domain sample can be explained quite well. This S/O 
pattern-dependent hopping model was successfully ex- 
tended to quantitatively explain the reflectivity spectra 
R(ui) of other NSMO multi-domain crystals. This clearly 
demonstrates that the polaron scenario is indeed valid, 
and that the polaron hopping motion should be strongly 
affected by the S/O ordering patterns. 

The NSMO single crystals were grown by the floating 
zone method [ fhi] . After polishing and proper annealing, 
R(lu) were measured by near normal incident unpolar- 
ized light over a wide photon energy region 5 meV— 30 
eV. For the single-domain a;=0.55 crystal, the crystal axes 
were determined using the neutron scattering measure- 
ment [Q. Using polarized light, we obtained R(u>) for 
the aft-plane, i.e. R a b(u>), and for the c-axis, i.e. R c (u>). 
Details of the optical measurements were described in 
Ref. g|. 

Figure 1(a) shows R(oj) for NSMO multi-domain sin- 
gle crystals at room T. Note that all the samples are 
insulators with disordered S/O states at room T |]l3| , |l6| . 
Clearly, R(u>) are similar for all the samples despite the 
different hole concentrations. On the other hand, Fig. 
1(b) shows R(u}) for the NSMO multi-domain samples 
at low T. R(uj) show quite different characteristics for 
the samples with different doping level x. For the x=0A0 
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(i.e., F-type ordered) sample, R(cu) increases significantly 
below 1.2 eV and approach 1.0 in the dc limit, show- 
ing metallic behavior For the :r=0.50 (i.e., C-E-type 
ordered) sample, R(u>) shows an insulator-like behavior 
with rather strong suppression below 0.2 eV. For the 
:e=0.55 (i.e., .A- type ordered) sample, R(ui) keeps increas- 
ing and approaches to about 0.6 as u> decreases. Finally, 
for the 2=0.65 (i.e., C-type ordered) sample, R(oj) shows 
an insulating behavior. Such large differences in R(uj) at 
low T are not easily understood by keeping in consider- 
ation of the similarity in R(u>) at room T and the small 
differences in x. It is quite apparent that S/O ordering 
patterns could be playing an important role in optical 
properties of NSMO. 

We have obtained c(w) using the Kramers-Kronig 
(KK) analysis of R(u>) at room T. Since the optical 
anisotropy in these disordered S/O states is expected to 
be small, the KK analysis can be applied at room T even 
for the multi-domain crystals. As shown in Fig. 2(a), 
o~(lj) are quite similar for all the samples. In order to 
compare the polaron and the orbital excitation scenar- 
ios, we should obtain terms related to the coherent and 
the incoherent motions of free carriers from o~(u>). 

In the perovskite manganites, there are numerous op- 
tical transitions. As shown in Fig. 2(b) for the :r=0.55 
sample, the room temperature cr(w) below 4.0 eV was fit- 
ted with four peaks represented by the dotted lines. The 
first peak (Peak I) was due to either polaron hopping 
motions [ pl]Jl2] | or orbital excitations |Q. The second 
peak (Peak II) located around 1.5 eV was suggested to 
be due to the transition between JT split e g bands Jl7| . 
[A Monte Carlo simulation using a Hamiltonian with the 
double exchange and JT interactions also clearly pre- 
dicted the existence of Peaks I and II |[2]].] The third 
peak (Peak III) was suggested to be due to a transition 
between the Hund's rule exchange ( J ~ 2 eV) split bands 
fill ■ Finally, the fourth peak (Peak IV) was proposed to 
be due to the charge transfer transition between the O 2p 
and Mn 3d levels . Note that the position of each peak 
are quite consistent with the results of other measure- 
ments jl9| , and the validityof those peaks has been sug- 
gested by many groups j|,[l5|,^0| and such analyses have 
already been performed quite successfully to explain the 
optical responses of (Bi,Ca)Mn0 3 @, Nd .7Sro.3Mn0 3 
|| , and (La,Ca)Mn0 3 @. 

In the polaron scenario, Peak I should be due to the JT 
polaron hopping motion. Since it can be considered as 
the interatomic transition from the JT distorted Mn 3+ [x 
ions per unit volume] to Mn 4+ [(l-i) ions per unit vol- 
ume], it should be proportional to x(l— x). The inset 
in Fig. 2(a) shows the x-dependent strength of Peak I 
(Si). The dotted line in the inset shows, within our fit- 
ting error bars, that the ^-dependence of Si is well fitted 
with x(l—x). A similar behavior was also observed in 
(La,Ca)Mn03 jo], supporting the agreement with the 
polaron scenario. However, as far as we know, the in- 



dependence of Peak I in the orbital scenario, i.e., the 
incoherent motion of free carriers due to orbital excita- 
tions, has not been predicted yet. 

To obtain further insights, we have investigated o~{ui) 
for the i=0.55 single-domain crystal at low T in the 
A-type S/O ordered state. The open triangles in Figs. 
3(a) and 3(b) show er(cj) along the c-axis, cr c (u)), and in 
the aft-plane, u a b(ui), respectively. Clearly, there exists 
a rather strong anisotropy of a(uS) below 1.5 eV. The 
solid lines show the fitting curves based on the above- 
mentioned mode assignments. As shown in the dotted 
lines, the strengths of Peaks II and IV were assumed to 
be T- independent, and the strength of Peak III were as- 
sumed to be nearly zero. Since Peaks II and IV are due 
to the interband transitions between JT split e g bands 
and the charge transfer transition between the O 2p and 
Mn 3d bands, respectively, they should be nearly T- 
independent. And, since Peak III is due to a transition 
between the Hund's rule exchange ( J ~ 2 eV) split bands 
[jUll , it should vanish in the spin ordered states. On the 
other hand, the strength of Peak I is suppressed along 
the c-axis but increases in the aft-plane. The enhance- 
ment factor of Peak I in the A- type S/O ordered state 
along the c-axis, rj A , was estimated to be about 0.4, and 
that in the aft-plane, r] A b , about 1.8. The difference be- 
tween rj A and r) A for Peak I can explain the reported op- 
tical anisotropy in the S/O ordered state. It is important 
to note that all the parameters, except the enhancement 
factors, are fixed in the fitting process: namely, the solid 
lines in Figs. 3(a) and 3(b) are the results of single pa- 
rameter fittings, keeping all other parameters the same 
as in the room temperature o~(oj) fitting. 

How can we explain the difference between rj A and rj A b ? 
Recently, <r(ui) for the planar orbital ordered phase of 
the e g electron with x 2 — y 2 symmetry (the A- type orbital 
ordered state) were calculated based on the orbital exci- 
tation scenario M] . This 2-dimensional model predicted 
that the spectral weight of the incoherent absorption, 
Si, should decrease with decreasing T. On the other 
hand, that of the coherent absorption, corresponding to 
the Drude peak, should increase. Contrary to the pre- 
dictions, the experimental value of rf A becomes larger 
than 1.0. Moreover, Fig. 3(b) does not show any signif- 
icant Drude-like absorption term. [Note that the value 
of a a b(uj) at the dc limit in Fig. 3(b) was estimated at 
around 2xl0 2 f2 -1 cm -1 , which is consistent with the 
transport experiment value on a single-domain crystal 
Jill]. Therefore, the model based on the orbital excita- 
tion scenario cannot explain our experimental data. 

As far as we know, there are no explicit predictions on 
how the S/O ordering pattern will affect the polaron ab- 
sorption. However, as shown in Fig. 3(c), it is quite clear 
that the polaron hopping motion along the c-axis should 
be suppressed since the antiparallel spin alignment in the 
A- type S/O ordered state will cause a large energy cost 
of J. Due to the x^ — y 2 e g orbital ordering within the 
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a&-plane, the overlap of two Mn e g orbitals via an oxy- 
gen p z orbital along the c-axis becomes small, thereby 
suppressing Si even further. On the contrary, Si in the 
a&-plane will become enhanced due to the ferromagnetic 
spin alignment and the increase of orbital overlap in the 
j? — y 2 orbital ordered state. 

In the spin/orbital pattern dependent polaron 
(SOPDP) model, the changes in Si in the S/O ordered 
state should be calculated by considering the orbital con- 
tributions in the conventional polaron models, which usu- 
ally include the double exchange and the JT polaron 
terms. However, this calculation might be quite difficult. 
In this paper, we made simple and crude estimations of 
the enhancement factors using the matrix element terms 
in the dipole transition. First, we considered the JT po- 
laron absorption as a transition from a lower e g level of 
a JT distorted Mn 3+ site to one of doubly degenerate 
e g levels of Mn 4+ sites. Second, for the ^4-type ordered 
state, we evaluated the sum of possible dipole matrix el- 
ements considering the neighboring S/O configurations. 
For the Mn 3c? e g and the O 2p orbital wave functions, 
we used the wave functions of isolated ions located at 
corresponding sites. Third, for the disordered state, we 
evaluated the sum of possible dipole matrix elements, as- 
suming that there is no preference of S/O for the initial 
state. Finally, we divided the summed values to obtain 
enhancement factors of Si in the ordered state. Details 
of this estimation procedure will be published elsewhere 
@. We found that rf A ~ and rjf ~ 1.7. The estimated 
value of r\ a A h agrees quite well with the experimental value, 
suggesting the validity of the SOPDP model. The sup- 
pression of Si along the c-axis can be explained at least 
qualitatively. However, the difference between the exper- 
imental and the estimated values of rf A might be due to 
the experimental errors and/or the over-simplification of 
our estimation. 

For the other multi-domain single crystal samples, the 
polarized optical microscope images showed that they 
consist of multi-domains whose typical size is a little 
larger than 10 fim. Thus, the KK analysis based on R(u>) 
for multi-domain crystals may lead to errors in the posi- 
tions and strengths of Peaks I and II ||l]] . To avoid such 
errors, we decided to analyze our R(lo) directly. The op- 
tical responses of inhomogeneous media have been inves- 
tigated extensively. In the long wavelength limit, where 
the wavelength A of the incoming light is larger than the 
domain size L, the medium should look homogeneous 
with an effective conductivity (22] . On the other hand, 
if A < L, the light will see the responses of individual 
domains, so we will obtain averaged optical responses. If 
the c-axes of the domains are randomly oriented along 
any one of the crystallographic axes, the measured R(u>) 
is given by [%R ab (w)+±R c (w)]. 

To test the above argument, we have calculated R(lj) 
for the multi-domain case of A-type (i.e., ir=0.55) or- 
dering, whose polarization-dependent <j(oj) were obtained 



experimentally. As mentioned above, the predictions of 
the SOPDP model are quite well matching to the real 
o~(uj) of a single-domain crystal. From a a b(uj) and <7 c (lu) 
of the SOPDP model, as shown in Figs. 3(a) and 3(b) 
(solid lines), respectively, corresponding reflectance spec- 
tra for given polarization directions were obtained. Using 
the predicted R a b{u) and R c (uj), R(w) were evaluated 
and is shown as the open triangles in Fig. 1(b). The 
excellent agreement between the experimental data and 
the SOPDP model predictions suggests that the method 
of obtaining the R(lu) for a multi-domain sample using 
the SOPDP model is accurate for A < L ~ 10 \xm. 

For the .F-type (i.e., i=0.40) sample, which has 
isotropic S/O structure even at low T, the enhancement 
of the mid-IR region and the metallic behavior can be ex- 
plained well with the existing JT polaron model, which 
considers the JT distortion and double exchange interac- 
tion, as observed in Lao.7Cao.3Mn03 0. For other multi- 
domain samples, similar analyses based on the SOPDP 
model were also performed |2l]| . For the C-type ordered 
(i.e., i=0.65) sample, rf c ~ 2.9 and v$ ~ can provide 
a good description to our experimental data. [Here we 
used the same method with t]a using 3^-r 2 orbital for 
the initial state instead of x 2 -?/ 2 at low T] For the CE- 
type (i.e., x=0.50) sample , this approach can provide a 
reasonable fit to the experimental R{lo) when Peak I be- 
comes asymmetric 21 . It might be due to the opening 



of the charge gap at the low frequency region in the S/O 
ordered state. As shown in Fig. 1(b), most experimen- 
tal R(u>) can be explained by the theoretical predictions 
using the SOPDP model with a single fitting parame- 
ter (77). The good agreements suggest that the optical 
response observed in various spin/orbital ordered states 
of NSMO can be explained quite well within the SOPDP 
model. 

In conclusion, we have presented the temperature- 
dependent spectral weight changes of Ndi-^Sr^MnOa 
single crystals, which strongly depend on the pattern of 
the spin/orbital ordering at low temperature. These op- 
tical responses can be described quantitatively by taking 
account of spin/orbital pattern-dependent polaron hop- 
ping. 
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FIG. 1. R(lj) of NSMO at (a) 300 K and (b) 15 K. In (b), 
solid circles, solid squares, open triangles, and open circles are 
the predictions of the SOPDP model for x=0.40, 0.50, 0.55, 
and 0.65, respectively. 
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FIG. 2. (a) a(tu) of NSMO at room T. The inset shows 
^-dependent Si. The solid circles represent data from this 
study, and the open circle is from the published data by Lee 
et al jlf| . The dotted line in the inset represents the functional 
form of x(l-x). (b) a(to) of NSMO (z=0.55) at room T (solid 
circles), where the solid line and dotted lines represent the 
fitting curves. 
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FIG. 3. The experimental (open triangle) and SOPDP 
model-predicted <j(lo) (solid line) of NSMO (x=0.55) at low 
T for (a) the afe-plane and (b) the c-axis. (c) Schematic rep- 
resentation of the SOPDP model in the A-type ordered state. 



